Prediction of new crystal structure phases in metal borides: a lithium 

monoboride analog to MgB2. 
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Modern compound prediction methods can efficiently screen large numbers of crystal structure 
phases and direct the experimental search for new materials. One of the most challenging problems 
in alloy theory is the identification of stable phases with a never seen prototype; such predictions 
do not always follow rational strategies. While performing ah initio data mining of intermetallic 
compounds we made an unexpected discovery: even in such a well-studied class of systems as metal 
borides there are previously unknown layered phases comparable in energy to the existing ones. 
With ah initio calculations we show that the new metal- sandwich (MS) lithium monoboride phases 
are marginally stable under ambient conditions but become favored over the known stoichiometric 
compounds under moderate pressures. The MS lithium monoboride exhibits electronic features 
similar to those in magnesium diboride and is expected to be a good superconductor. 



Development of theoretical methods able to guide the 
experimental search for new materials with desired prop- 
erties is a fast-growing field of research in materials 
sciencelH |lllllll|6[[lli. A particular effort has 
been put into the coupling of ah initio electronic structure 
methods with efficient data-mining algorithms to deter- 
mine and utilize correlations in binding mechanisms in a 
chosen set of structures 0, 0, 0; 0- However, even with 
the most advanced optimization algorithms one has to 
restrict the search space to a given structure 4], lattice 
typej^, or prototype libraryjsT Q|- Moreover, in some 
approaches structures are kept symmetry-constrained or 
not relaxed at all. Considering the rich nature of bond- 
ing in solids, such limitations can lead to overlooking the 
most stable phases. This makes identification of new pro- 
totypes a vital step towards a more complete description 
of alloys. 

Expanding the library of ah initio energies of binary 
alloys described in Ref. we observed that one of the 
phases in the Mg-B system, A2B2 fcc-(lll) (or V2 0), 
unexpectedly underwent significant structural relaxation 
and became comparable in energy to the mixture of two 
coexisting stable phases, MgB2 and hcp-Mg. Upon ex- 
amination of the relaxation process we found that there 
is a continuous symmetry-conserving path from V2 (Fig. 
Ic) to a new metal-sandwich (MS) structure MSI (Fig. 
lb) The latter has four atoms per unit cell with 
hexagonal layers of boron separated by two triangular 
layers of metal; it bears strong resemblance to the AIB2 
prototype (Fig. la) having a similar sp^ boron-boron 
bonding. The extra metal layer is inserted in a close- 
packed fashion with Mg-Mg bond length close to that in 
the pure hep structure, which may account for the near- 
stability (by a few meV/atom) of the magnesium mono- 
boride. To the best of our knowledge, this structure has 
not been considered before for any binary alloy. We have 
constructed a series of related MS structures with differ- 
ent stacking sequences of metal and boron layers (such 
as MS2 in Fig. Idlll), which is a generalization of the 



MSI prototype. 

The identification of the MS structures has prompted 
us to re-visit metal borides, which have lately been a 
focal point in the search for new superconducting ma- 
terials, following the discovery of superconductivity in 
MgB2 with a surprisingly high transition temperature 
of 39 K 12J. The high T^ in MgBa was explained as 
due to a very strong electron-phonon coupling of the E2g 
mode to the hole-like pcr-states of boron 13] . These states 
are also crucial for structural stability of metal diborides 
and their availability in these compounds implies weak 
binding^3|- This feature makes the effectively electron- 
doped (because of the extra metal layer) MSI structure 
a particularly good candidate to be a stable configura- 
tion for low-valent metal borides. We have investigated 
a large set of alkali, alkaline and transition metal borides 




a) AIB2 
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FIG. 1: (color online). Known (AIB2, V2) and proposed 
metal-sandwich (MSI, MS2) prototypes. The hexagonal lay- 
ers of boron (grey) are separated by triangular layers of metal 
(yellow, only four atoms per layer are shown). 
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with ab initio calculations and found that MSI and MS2 
phases in one system, Li-B, are indeed stable enough to 
compete against known phases. 

Present ab initio simulations have been performed 
with VASP nil H^ : we use projector augmented 
waves (PAW) [ig and exchange-correlation function- 
als as parametrized by Perdew, Burke, and Ernzerhof 
(PBE)'ls| for the Generalized Gradient Approximation 
(GGA). For comparison purposes we carry out tests in 
the Local Density Approximation (LDA)[l9(, and with 
ultrasoft (US) pseudo pot entials(2Qi| . We also employ an 
LAPW code WIEN2K to plot characters of electronic 
bands and partial density of states (PDOS). All struc- 
tures were fully relaxed. Numerical convergence to within 
about 2 meV/atom is insured by high energy cut-off (398 
eV) and dense k-meshesj2^. 

Despite the potentially important application to bat- 
teries, the Li-B system is not fully explored and the 
question of the most stable monoboride phase remains 
open|3 EE HE H^' Successful experimental attempts 
to synthesize lithium monoboride under ambient pres- 
sure date back to the 1970s j2^ (29j. However, the struc- 
ture of g-LiB, shown in Fig. 2a, has been identified only 
recently HE Hj, HEl- ^ theoretical study suggests that 
lithium atoms in the lithium monoboride could be sub- 
ject to low-energy motion along the linear boron chains, 
as a-LiB and a related /3-LiB phase (NiAs prototype, see 
Feg. 2b) were found to be nearly degenerate |2^. 

We compare our MS lithium monoboride to the known 
stoichiometric phases directly, using the experimentally 
observed a-LiB as a reference point. The following re- 
sults may contain relatively big systematic errors (i.e. 
due to approximations in the exchange-correlation func- 
tional and in the treatment of the charge density), since 
we consider structures with different types of boron- 
boron bonding. Wc find that the MSl-LiB phase is 
3 meV/atom lower in energy than a-LiB and about 7 
me V/ atom higher than /3-LiB[23- The a-LiB and /3- 
LiB phases differ by a noticeable 10 meV/atom, a sur- 
prising result considering that our fully-relaxed cell pa- 
rameters (a-LiB: a = 4.026 A, c = 3.112 A and /3-LiB: 
a = 4.010 A, c = 3.119 A) closely match the reported 
LDA values |2a|. We check the sensitivity of these results 
to the approximation used by carrying out two more tests 
with WIEN2K (LAPW, in the GGA-PBE) and VASP 
(PAW-LDA). 

The relaxed cell parameters in the LAPW-GGA calcu- 
lations are nearly the same as in the PAW-GGA 10, llj. 
In this approach MSl-LiB is a little lower in energy 
(13 meV/atom below a-LiB) and the most stable of the 
three 32] . The PAW-LDA gives a much lower relative en- 
ergy for MSl-LiB (42 meV/atom below a-LiB). A big dis- 
crepancy in the description of relative energies with the 
two approximations has also been seen in carbon systems 
with different types of bondingjs^; the GGA is expected 
to be more accurate. In all the tests above a-LiB and 
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FIG. 2: (color online). Calculated relative enthalpy, 
AH(P)=H(P)-Ha_z,iB(P), as a function of pressure for 
lithium monoboride phases. The points are connected with 
splines. The larger yellow spheres are lithium atoms and the 
smaller black spheres connected with sticks are boron atoms. 



/3-LiB remain 10 meV/atom apart. The independence of 
this result on the calculation approach can be explained 
by the similarities of the two structures, a trend also ob- 
served in other systems !33|. 

We conclude that within the accuracy of these calcu- 
lations the MSl-LiB phase is comparable in energy to 
the known phases. Synthesis conditions could therefore 
be the deciding factor for which compound forms in the 
Li-B system at 50% concentration. We next evaluate the 
effect of hydrostatic pressure on the relative stability of 
the lithium monoboride phases. The three phases have 
very close volumes per atom at zero pressure. However, 
the a-LiB and /?-LiB with linear boron chains are much 
harder than the sp'^-bound MSl-LiB: the respective vol- 
umes at P = 12 GPa are about 85%, 85%, and 69% 
of the zero-pressure values (interestingly, the intralayer 
bonds in MSl-LiB undergo non- monotonous expansion 
with pressure, stretching by about 2% at 5 GPa). The 
big compressibility of the MSl-LiB is an expected feature 
for a layered structure with a weak interlayer binding: 
at P = 12 GPa the c-axis shrinks down to 68% of its 
zero-pressure value, primarily due to the decrease in the 
Li-Li interplanar distance. The weakening of the metal- 
lic binding in the lithium bilayer is apparently caused 
by the significant charge transfer from lithium to the 
boron layer. The extreme softness of MSl-LiB suggests 
use of high pressures for synthesis of this phase. We 
plot the calculated enthalpies of the considered lithium 
monoboride phases with respect to a-LiB versus external 
pressure m Fig. El The plot shows that even if the MS- 
LiB are not the most stable phase at ambient pressure, 
they are quickly driven below the known phases and be- 
come favored by over 90 me V/ atom at P 12 GPajs^. 
Note, that the hydrostatic pressure removes the degen- 
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FIG. 3: (color online). Slices of the potential energy sur- 
face in MS2-LiB along the eigenvectors of the softest phonon 
modes at V . The modes represent the interlayer sliding in the 
X (see the [010] view) and y directions of undistorted boron 
layers rigidly locked with the adjacent lithium layers. The in- 
set demonstrates two shallow minima at small displacements. 
Ah initio data points are connected with splines. 



eracy of the MSl-LiB and MS2-LiB, different only by a 
long-period stacking shift: the pressure favors MS2-LiB, 
in which the flatness of the boron layers is enforced by 
symmetry. Use of appropriate substrates could also pro- 
mote growth of the new layered structures. 

An important question is how to detect the new phases 
if they indeed form. Because the MSI and MS2 are nearly 
degenerate at low pressures, the resulting structure at 
50% concentration in Li-B could be a random mixture 
of the two and would be hard to detect with standard 
x-ray methods. Fortunately, the two phases have similar 
interlayer distances for boron and one would expect to 
see an x-ray peak (at A — 1. 5418 A, zero pressure) in 
the range from 2Q = 16.1° for MS2 to 26* = 16.6° for 
MSI . None of the observed peaks in the samples prepared 
at ambient pressure |2^ match these calculated values. 
Experimental effort to synthesize lithium monoboride at 
high pressures is highly desired. 

Ah initio studies have shown that existence of ener- 
getically favorable phases could be barred by dynami- 
cal instabilities 36] . We have investigated phonon modes 
at F-point with a force-constant matrix diagonalization 
method, implemented in VASP, using the MS2 structure 
with a convenient rhombohcdral unit cell. MS2-LiB has 
two weakly interacting boron layers per unit cell which 
leads to an effective double-degeneracy of the known 
phonon modes at F in the AlB2-type compounds. Three 
of the remaining six optical modes, which involve sliding 
of undistorted layers, are rather soft. The soft modes are 
defined by the weakened interlayer Li-Li interaction and 
the perfectly stacked structure may actually be prone to 
small interlayer shifts, as there are two shallow minima 



at about 2% displacements of the in-plane lattice vector 
a (see Fig. 3). We crudely estimate the frequency of 
the soft modes at F by using all points up to 5% dis- 
placements and obtain ijjx^y=hh cm~^ and Wz=86 cm"-'^. 
The modes at the A point are slightly softer [u}x,y—'il 
cm^^ and ujy^—h'd cm^^) but they remain real (interest- 
ingly, the frequencies are comparable to those in another 
layered structure, graphite |33). The interlayer sliding 
phonon modes along the F-A direction are expected to 
be the softest optical phonon modes in the whole Bril- 
louin zone: any finite k^^y would result in distortions of 
the rigid boron or lithium layers and a consequent phonon 
hardening. These considerations indicate that the MS2- 
LiB phase is dynamically stable. Hydrostatic pressure 
is again found to be a stabilizing factor: already at P=2 
GPa the small distortions are suppressed and the frequen- 
cies of the softest x-y phonon modes are nearly doubled. 
Calculation of phonon modes for the full Brillouin zone is 
a subject of future study. With all the phonon branches 
calculated one can also estimate the vibrational entropy 
contribution to the Gibbs free energy to complete the 
analysis of the thermodynamic stability of MS-LiB. 

The investigation of the phonon modes reveals that the 
actual stable unit of the MS-LiB phases is a hexagonal 
boron layer coated with two lithium layers on each side 
and can be related to graphene. Indeed, according to 
our estimates the cohesion per area between these B2Li2 
sheets and their curvature energy are comparable to those 
for a carbon layer [s^. One could speculate that if MS- 
LiB was exfoliated into such sheets, they might form 
carbon-like tubular or nanoporous structures. It would 
be interesting to see whether stable MS compounds could 
be obtained in other intermetallics based, for instance, on 
beryllium. Going beyond binaries also seems a promis- 
ing direction to find more stable compounds: additional 
metal sites give more flexibility to experiment with dif- 
ferent ternary alloys. 

Finally, we give a brief discussion of the electronic 
properties of the MS phases focusing primarily on the p- 
states in boron, important for superconductivity in metal 
diborides. Our LAPW calculations [U |^ show that 
the boron per band along F-A in MS2-LiB has practi- 
cally no dispersion (Fig 3), an expected feature caused 
by the large separation between boron layers. Com- 
pared to a hypothetical AlB2-LiB2 phase the F-A band 
downshifts by 0.6 eV and the ipa PDOS at Fermi level 
drops by 25% to iV^^(O) = 0.046 states/ (eV-spin) per 
boron atom, which indicates higher flUing of these bond- 
ing boron states. Note, that bonding p7r states also be- 
come occupied (Ej? is located exactly in the middle of 
the pTT pseudo-gap). The combination of these effects 
explains the dramatic stabilization of the lithium mono- 
boride. Remarkably, the pa PDOS in MS2-LiB stays 
about 12% higher than that in MgB2. 

We flnd the frequency of the anharmonic E2g mode 
at F-point (649 cm"-'^) in MS2-LiB to be quite close to 



4 




r M K TA L H r 0.1 0.0 

k space PDOS 



FIG. 4: (color online). Band structure and partial de nsity 
of states (PDOS) in MS2-LiB, calculated in LAPWlU E^ . 
PDOS units are states/(eV-spin) per boron atom. The thick- 
ness of band structure lines is proportional to boron px,y (red) 
and pz (green) characters. 



that (624 cm""'^) in MgB2[23j which can be related in 
part to the almost perfect match of the covalent boron- 
boron bond lengths The similarities of the boron 
layer properties in these borides appear to extend to the 
electron-phonon coupling as well: following procedure in 
Ref. |42|, we estimate the deformation potential V for 
the Fi2g mode to be 13 eV/A in both phases. Unlike the 
two-band superconductor MgB2, MS2-LiB has negligible 
contribution from the boron pir states at Ep- These sim- 
ple arguments suggest that the superconducting Tc in the 
MS2-LiB would be at least comparable to that in MgB2. 

In conclusion, this work demonstrates that our knowl- 
edge of binary alloys is still incomplete and the identi- 
fication of yet unknown phases may require new predic- 
tion strategies. The rather accidental finding of the MS 
phases should be credited to the exhaustive considera- 
tion of all candidates used in our data-mining approach 
and the careful structural relaxation in the calculation of 
their ground states. The MS structures are shown to be 
particularly suitable for the Li-B system: lithium mono- 
borides are marginally stable under ambient conditions 
but become favored over the known stoichiometric com- 
pounds under pressure. The effectively electron-doped 
MS-LiB phases retain a significant PDOS from the pa- 
boron states a,t Ep and could be the long-sought non- 
magnesium-based high Tc metal boride. Moreover, rela- 
tive to MgB2 the MS-LiB phases are demonstrated to be 
hole-doped; if synthesized, they would provide valuable 
information on how the hole-doping affects the Tc in sto- 
ichiometric metal borides and whether MgB^ really has 
the highest possible Tc in its class liilil lHlifij i. 

We thank F.H. Cocks, V. Crespi, P. Lammert, R. 
Margine, and J. Sofo for valuable discussions. 
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